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We construct a novel observable for finite temperature QCD that relates confinement and chi- 
ral symmetry. It uses phases as boundary conditions for the fermions. We discuss numerical 
and analytical aspects of this observable, like its spectral behavior below and above the critical 
temperature, as well as the connection to chiral condensate and center symmetry. 
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The main phenomena in QCD at the finite temperature transition are deconfinement and chiral 
symmetry restoration. There has always been the question whether there is a mechanism connec- 
ting the two, in particular since in the quenched theory the corresponding phase transitions occur 
at the same critical temperatures [1]. Here we define an observable that indeed links the two and 
discuss numerical findings from quenched lattice configurations [2] . 

The Polyakov loop, on the one hand, is the order parameter of confinement, being traceless in 
the confined phase and moving towards the center of the gauge group at temperatures above Tc. The 
spectral density of the Dirac operator at the origin, p (0) , on the other hand, is the order parameter 
of chiral symmetry. It is related to the condensate by the famous Banks-Casher relation [3] {'^ff^f) = 
— 7rp(0) and vanishes above Tc. How does confinement leave a trace in the Dirac spectrum? The 
answer lies in the dependence on temporal boundary conditions, as we will show now. 

We use the lattice as a regulator. The (untraced) Polyakov loop is the product of temporal 
links. For the lattice Dirac operator we use the staggered one [4], which can be viewed as hopping 
by one link. It is well-known that the fc-th power of the Dirac operator at the same argument, 
D^{x,x), contains all products of links along closed loops of length k starting at x. To distinguish 
the Polyakov loop from 'trivially closed' loops (like the plaquette), one needs phase boundary 
conditions for the fermions [5], ^{xq + j3,f) = e"l'\i/{xo,x). These boundary conditions amount to 
an imaginary chemical potential. They can be easily implemented by replacing the temporal links 
Uq in some time slice by e'^Uo. As a consequence, all loops get a factor e'l"' according to their 
winding number q around the temporal direction, trivially closed loops remain unchanged. 

In this way, the Polyakov loop can be reconstructed from the Dirac spectrum by using at least 
three boundary conditions, see [6]. This 'thin' Polyakov loop, however, has poor renormalization 
and scaling properties and it turned out that in this approach it is UV dominated [6]. 

Influenced by the Jena group [7], we instead consider the propagator l/(wi + D) with some 
probe mass m and use a geometric series obtaining all powers of the Dirac operator. Denoting the 
Dirac operator at a particular boundary condition <p by D(p, i.e., including the factors of the 
propagator is given as a product of links along all closed loops I, 

where |/| is the length of the loop and sign(/) comes from the staggered factor. 

Of importance in ( 1 ) is the phase factor, since one can project onto a particular winding ^ by a 

Fourier transform w.r.t. the boundary angle, J^'^ dcpe^'f"^ Specifying to a single winding, q = l, 

like for the Polyakov loop, we arrive at [2] 

r-^ / _^ ^(tr.n W). (2) 

Jo 2k v\ m + D^I mV^(^j(2am)l'l \ ( ' / 

We refer to the new observable £ as the 'dual condensate', because it is obtained thi^ough a Fourier 
transform from the trace of the propagator. Indeed, in the massless limit (after the infinite volume 
limit as usual) we obtain the chiral condensate integrated with a phase factor over the boundary 
conditions and, using the Banks-Casher relation at every individual angle (p, the corresponding 
representation in terms of the eigenvalue densities p(0),j,. The right hand side of (2) represents the 
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Figure 1: Left: The dressed Polyakov loop E as a function of T at various lattices (for m — 100 MeV). 
Right: The corresponding plot for the conventional Polyakov loop. 



'dressed Polyakov loop', that is the set of all loops which wind once around the temporal direction. 
In the infinite mass limit, detours become suppressed and only the straight Polyakov loop survives 
as it is the shortest possible loop in this set. 

Fig. 1 shows that £ is indeed an order parameter. Keeping the mass m fixed, £ vanishes 
below the critical temperature (which is about 280 MeV in the quenched case) and develops an 
expectation value for higher temperatures. One finds that the results (when expressed in physical 
units) obtained for different volumes and with different resolution essentially fall on a universal 
curve. This illusti'ates the good renormalization properties of our observable, which are inherited 
from {i^f^) and may be viewed as an effect of the dressing rendering such loops less UV dominated. 

In a spectral representation, the Dirac operator in E is replaced by an eigenvalue sum. 

As the eigenvalues appear in the denominator, we expect the sum to be dominated by the IR modes. 
As Fig. 2 shows, this is confirmed by our lattice data. 
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Figure 2: Individual and accumulated contributions to the spectral sum (3) at two different masses. 
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Figure 3: The integrand (I,(m + 4'V>/^ 

as a function of the boundary condition ^ for two values of am 
and two temperatures, for configurations with real Polyakov loop. 

How is a finite resp. vanisliing order parameter £ built up by the eigenvalues? Fig. 3 shows 
that they respond differently to the boundary conditions in the confined vs. deconfined phase. The 
eigenvalues are independent of the boundary condition in the confined phase, which leads to a 
vanishing order parameter £. In the deconfined phase, on the other hand, the eigenvalues show a 
typical cosine-type of modulation. Together with the Fourier factor this yields a nonvanishing £. 

The chiral condensate has to behave essentially in the same way, as it is the integrand in the 
massless limit. Although the chiral condensate is finite in the confined phase, it is independent 
of the boundary condition <p and hence results in a vanishing dual condensate. The situation in 
the deconfined phase might be more confusing at first glance as the spectrum has a gap there. For 
boundary conditions in Une with the original Polyakov loop, however, the chiral condensate persists 
above the critical temperature T^, which has been seen in several studies [8]. This mechanism is 
needed to ensure a finite £ and should actually be at work for all T > Tc. 

In the quenched case the conventional Polyakov loop is the order parameter for center symme- 
try. Under center transformations the dressed Polyakov loop behaves in the same way. Hence, L is 
an order parameter for center symmetry, which is underneath our numerical findings and which we 
therefore expect to be independent of the choice of the lattice Dirac operator. 
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